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Abstract 


Jraphs  of  expected  dosage  from  an  infinite  line  source  are  presented  for  various 
wind  speeds,  heights  of  release,  and  stability  conditions.  The  graphs  are  based  on 
Sutton's  diffusion  equation  with  parameters  derived  from  Project  Prairie  Grass  data 
by  Haugen,  Barad,  and  Antanaitis..  A  method  of  estimating  the  range  of  validity  of 
these  equations  for  a  line  source  of  finite  length  is  also  presented. 
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Operational  Prediction  of  Diffusion  Downwind 

From  Line  Sources 


1.  INTRODUCTION 

The  Air  Force  has  a  requirement  for  knowledge  of  the  ground-level  distribution 
of  various  gases  and  aerosols,  released  in  a  quasi-instantaneous  fashion,  along  lines 
that  may  either  be  ground-based  or  elevated.  The  purpose  of  this  Survey  is  to  pro¬ 
vide  a  method  whereby  operations  that  require  such  information  can  be  planned . 
Graphs  of  the  expected  ground-level  dosage  divided  by  source  strength  are  given  as  a 
function  of  distance  downwind  from  the  release  line  for  various  heights  of  release, 
wind  speeds,  and  the  stability  of  the  lower  atmosphere. 

The  equations  from  which  the  graphs  were  constructed  are  based  on  those  derived 
by  Sutton.*  While  other  equations  have  been  presented,  Sutton's  are  the  best  known 
and  most  widely  studied  and  require  fewer  specialized  measurements.  In  order  to 
appreciate  the  limitations  of  the  results,  presented  in  Appendix  A,  a  brief  review  of 
Sutton's  equations  is  in  order.  No  one  should  attempt  to  use  the  results  without  first 
reading  the  following  sections. 


2.  THE  EQUATIONS 

Sutton  derived  his  equations  from  statistical  considerations  of  turbulent  diffusion 
(Received  for  publication  20  February  1964) 
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from  an  instantaneous  point  source.  By  integrating  with  respect  to  time,  he  obtained 
an  equation  for  a  continuous  ground-level  point  source  and  this  latter  equation  is  the 
basis  for  generalizations  to  other  source  types. 

The  equation  takes  the  general  form 


X  =  1 

"5  ”  iro  d  f 
y  Z 


exp 


2  21 

-V 

2o‘  2o‘ 

-  y  zj 


(1) 


where  X  is  the  concentration  at  any  point  (x,  y,  z),  Q  is  the  source  strength,  a 
and  oz  are  the  standard  deviations  of  the  distances  traveled  by  a  particle  in  the 
crosswind  (y)  and  vertical  (z)  directions,  respectively.  The  downwind  direction (x) 
is  the  direction  in  which  the  mean  wind  (u)  is  blowing  and  does  not  appear  explicitly. 
The  values  of  and  are,  however,  both  functions  of  x,  reflecting  the  fact  that 
the  cloud  spreads  both  laterally  and  vertically  as  it  travels  away  from  the  source. 
This  equation  is  presumed  valid  for  sources  based  at  the  ground  in  an  atmosphere 
where  the  wind  does  not  vary  with  height.  Furthermore,  it  assumes  that  any  tracer 
material  which  hits  the  surface  "bounces  "  back  into  the  air,  that  is,  there  is  no 
deposition  on  the  ground  or  other  surface.  It  also  assumes  the  crosswind  and  vertical 
distributions  are  Gaussian  and  that  the  downwind  speed  of  translation  (u)  is  much 
greater  than  the  downwind  diffusion. 

The  problem  then  is  to  express  and  oz  as  functions  of  x  and  of  atmospheric 
stability.  Sutton  presents  some  theoretical  expressions  of  the  form 


2o2  =  C2  x2‘n 


2o2  =C2x2'n 
z  z 


(2) 


where  the  C's  are  functions  of  the  gustiness  and  stability  but  not  of  x,  and  n  is  ' 

presumed  related  to  the  wind  profile,  which  is  in  turn  also  a  function  of  stability. 

Equations  (2)  state  that  the  downwind  variation  of  ff,  as  represented  by  the 

term  x2  n,  is  the  same  for  both  the  vertical  and  crosswind  terms.  However,  Barad 
2 

and  Haugen  have  presented  an  analysis  of  some  data  for  continuous  point  sources 

q 

(Project  Prairie  Grass  )  which  indicated  that  different  values  of  n  (referred  to 
hereafter  as  n^  and  nz)  are  more  appropriate.  This  concept  was  used  in  constructing 
the  graphs  presented  here. 

Since  the  problem  of  direct  concern  is  an  instantaneous  line  source  rather  than 
a  continuous  point  source,  Eq.  (1)  has  to  be  transformed.  Integration  of  Eq.  (1)  in 
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the  croeswind  plane  from  y  *  -  «  to  y  *  +»  gives  the  appropriate  solution  for  a  con¬ 
tinuous  line  source.  Then,  if  the  concentration  X  is  replaced  by  the  total  dosage  D 
(units  of  quantity  times  time  per  unit  volume),  and  the  source  strength  Q  is  taken  as 
quantity  per  unit  length  of  release,  a  solution  for  an  instantaneous  line  source  is 
achieved.  Finally,  by  referring  to  the  method  of  images^  the  equation  can  be  made 
appropriate  for  a  source  whose  height  above  ground  is  h.  The  final  equation  is  then 


D 

5" 


2 


_  _<2-nz)/2 

V»CZ  u* 


exp 


(3) 


This  equation  then  gives  the  value  of  D/Q,  at  ground  level,  as  a  function  of  the 
source  height  h,  the  wind  speed  u,  and  the  downwind  distance  x,  provided  the 
parameters  Cz  and  nz  can  be  specified. 


3.  ESTIMATION  OF  PARAMETERS 

Haugen,  Barad,  and  Antanaitis4  have  presented  values  of  the  diffusion  par¬ 
ameters  Cy  Cz,  ny,  and  nz<  based  upon  a  careful  analysis  of  the  Prairie  Grass 
data.  As  stated  above,  these  parameters  are  functions  of  stability.  Haugen  et  al. . 
chose  the  Stability  Ratio  (SR)  as  the  quantitative  measure  of  stability.  This  par¬ 
ameter  is  defined  as 


SR 


X  105  °C  sec2 


- 2 


(4) 


The  subscripts  refer  to  the  height  above  ground,  in  meters,  at  which  the  tem¬ 
perature  (T)  and  wind  speed  (u)  were  measured.  This  quantity  is  negative  in  lapse 
situations  and  positive  in  inversion  situations. 

The  appropriate  values  of  the  parameter  for  various  stability  classes  as  found 
by  Haugen  et  al,  are  given  in  Table  1.  (The  integration  of  Eq.  (1)  with  respect  to 
y  eliminated  Cy  and  ny.  These  values  will  be  needed  later  in  this  paper,  however, 
and  are  presented  here  for  convenience.) 

The  values  appearing  in  Table  1  are  those  used  in  constructing  the  graphs. 
Obviously,  to  use  these  graphs,  one  must  first  be  able  to  ascertain  the  proper  sta¬ 
bility  category.  When  measurements  of  temperature  difference  and  wind  speed  at 
the  appropriate  levels  for  determining  the  stability  ratio  are  not  available,  the  follow¬ 
ing  may  be  useful  as  a  rough  guide  in  estimating  the  proper  stability  class. 

Based  on  the  experience  at  Project  Prairie  Grass,  during  daylight  hours  one 
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TABLE  1. 


Cs 

nz 

s 

n 

y 

Stability  Class 

Range  of  SR 

n  J2 

(meters) 

n  /  2 

(meters)  •' 

very  unstable 

<-0.45 

0.002 

-1.20 

0.38 

0.20 

moderately  unstable 

-0.45  to  -0.20 

0.02 

-0.40 

0.38 

0.30 

neutral 

0.20  to  0.20 

0.07 

0. 10 

0.38 

0.50 

moderately  stable 

0.20  to  0.45 

0.07 

0.20 

0.38 

0.65 

very  stable 

>  0.45 

0.07 

0.30 

0.38 

0.80 

would  expect  to  find  very  unstable  conditions  with  clear  skies  and  surface  winds 
less  than  10  knots;  neutral  conditions  would  be  found  with  overcast  skies  and  wind 
speeds  in  excess  of  15  knots.  If  neither  of  these  conditions  is  met  one  would  find 
moderately  unstable  conditions  during  the  day.  At  night  one  would  expect  very  stable 
conditions  with  clear  skies  and  wind  speeds  below  7  knots;  neutral  conditions  with 
overcast  skies  and  wind  speeds  greater  than  15  knots;  otherwise  one  would  expect 
moderately  stable  conditions.  In  addition,  near  Bunrise  and  sunset  one  expects  a 
period  of  neutral  conditions,  regardless  of  wind  and  cloud  conditions.  It  should  be 
emphasized  that  the  terms  "very*  and  "moderately"  are  based  on  the  conditions  dur¬ 
ing  Prairie  Grass.  Very  few  trials  were  made  with  surface  winds  less  than  3  knots. 
Therefore  one  might  find  that,  with  very  low  wind  speedB,  the  graphs  for  very  un¬ 
stable  and  very  stable  are  not  representative. 

A  question  concerning  the  wind  speed  to  use  in  entering  the  graphs  naturally 
arises.  It  is  well  known  that  the  wind  is  not  constant  with  height  as  the  theoretical 
development  presupposes.  Probably,  the  best  wind  speed  to  use  would  be  an  average 
between  the  surface  wind  speed  and  the  speed  at  the  release  height. 


4.  ACCURACY  OF  EQUATIONS 

It  must  be  stressed  that  values  of  the  parameters  used  in  this  study  have  been 
obtained  by  a  double  averaging.  They  were  derived  from  continuous  source  data 
which  is,  in  effect,  an  average  over  a  large  number  of  instantaneous  sources.  Also, 
the  values  of  die  parameters  represent  median  values  derived  from  a  number  of 
separate  continuous  source  emissions.  This  means  that  Eq.  (3),  with  these  par¬ 
ameters,  represents  what  would  be  expected  as  an  average  over  a  number  of  instan¬ 
taneous  releases  in  the  same  stability  conditions.  Equation  (3)  cannot  be  presumed 
to  be  more  than  an  educated  guess  for  any  particular  release. 
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No  data  from  elevated  sources  are  available  in  sufficient  quantity  or  quality  to 
determine  the  accuracy  of  the  graphs  in  Appendix  A;  only  subjective  estimates  of 
their  accuracy  can  be  made .  One  must  keep  in  mind  that  the  parameters  were  ob¬ 
tained  from  data  in  a  particularly  simple  geographic  location.  Furthermore,  the 
graphs  are  based  on  an  assumed  Gaussian  distribution  of  crosswind  and  vertical  con¬ 
centration.  The  assumption  of  a  Gaussian  distribution  seems  to  be  a  good  way  of 

r  ft 

representing  the  crosswind  distribution  from  a  point  source  but  there  is  evidence  ' 
that  the  vertical  distribution  may  be  better  represented  by  some  other  form.  How¬ 
ever,  neither  of  these  references  give  sufficient  information  from  which  the  appro¬ 
priate  diffusion  parameters  may  be  calculated.  Until  more  work  is  available,  it  is 
best  to  use  the  Gaussian  distribution. 

In  particular  applications,  one  would  expect  the  graphs  to  be  most  accurate  near 
the  distance  of  the  peak  dosage  and  when  the  source  is  close  to  the  ground.  At  dis¬ 
tances  greater  than  a  few  miles  from  the  source,  one  would  expect  deterioration  of 
the  results  because,  at  these  distances,  scales  of  atmospheric  motions  larger  than 
those  considered  in  ordinary  diffusion  equations  may  become  important.  Between 
the  elevated  source  and  the  maximum  dosage  the  lines  of  D/Q  vs.  distance  are  so 
steep  that  a  slight  shift  in  the  lines  could  result  in  a  difference  of  D/Q  of  several  or¬ 
ders  of  magnitude. 

In  nummary,  their,  one  might  expect  that,  around  the  maximum  value  of  D/Q, 
the  graphs  would  provide  estimates  of  D/Q  to  within  a  factor  of  4  about  2/3  of  the 
time  and  within  an  order  of  magnitude  almost  all  the  time.  The  estimated  distance 
from  the  source  to  the  maximum  D/Q  should  be  accurate  to  within  a  factor  of  2.  These 
statements  are  based  on  subjective  judgments  and  assume  an  accurate  determination 
of  the  source  height,  the  wind  speed,  and  the  stability  class. 

Another  point  which  must  be  stressed  is  the  relation  between  the  wind  direction 
and  the  direction  of  the  release  line.  The  derivation  of  Eq.  (3)  requires  that  the  wind 
direction  must  be  at  right  angles  to  the  release  line.  In  practice  it  appears  the  results 

would  not  differ  sufficiently  from  those  predicted  here  if  the  wind  direction  were  be- 

o  o  7 

tween  80  and  100  to  the  release  line.  Barad,  Haugen,  and  Fuquay  have  presented 

the  equation  necessary  to  consider  all  angles  between  wind  and  release  line  as  well  as 

finite  lengths  of  the  release  line.  However,  this  equation  can,  in  general,  only  be 

solved  by  digital  methods.  Solutions  of  this  equation  would  be  quite  valuable  and  should 

be  undertaken. 


S.  THE  LENGTH  OF  THE  RELEASE  LINE 

The  question  of  how  long  a  release  line  must  be  to  be  considered  effectively 
infinite  is  not  easy  to  determine.  The  farther  away  from  any  finite  source  the  more 
that  source  appears  as  a  point.  Conversely,  the  closer  one  is  to  the  center  of  a  finite 
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line  source  the  more  it  appears  as  an  infinite  line  source.  Thus,  one  would  expect 
some  region  down  wind  of  a  finite  line  source  to  be  unaffected  by  the  finite  length  of 
the  line.  This  concept  is  illustrated  schematically  in  Figure  1.  The  hatched  area  is 
the  region  down  wind  of  release  line  AB  in  which  the  finite  length  of  AB  is  immaterial 
for  diffusion  considerations.  The  problem  is  to  determine  the  boundary  of  this  region. 


Figure  1 

Consider  lines  AC  and  BD  normal  to  AB  and  line  CD  parallel  to  AB.  A  reasonable 
and  probably  conservative  assumption  can  be  made  which  should  provide  sufficient 
information  for  planning  purposes.  This  assumption  is  that  CE  and  FD  are  equal  to 
twice  the  crosswind  standard  deviation  ory.  Then  if  the  length  of  the  release  line  AB 
is  taken  as  L  and  the  length  of  the  line  EF  is  taken  as  L0>  it  follows  from  the  geom¬ 
etry  of  Figure  2  that 

TT  l2  =  20 y  •  (5) 

Substituting  for  cy  from  Eq.  (2) 

L  -  Lq  =  23/2  Cy  x(2'ny>/2  .  (6) 

The  value  of  Cy  is  not  dependent  on  stability  (see  Table  1)  and  is  equal  to  0. 38 
(meters)11^3.  Thus  with  sufficient  accuracy  Eq.  (8)  may  be  written,  if  L,  L0, 
and  x  are  measured  in  meters. 


L  -  L0  -  x1  ** 


(7) 
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where  the  ny  values  are  given  in  Table  1.  Therefore,  if  two  of  the  three  lengths  in 
Eq,  (7)  can  be  prescribed,  the  third  can  be  found.  Figure  2  shows  graphs  of  this 
equation  for  values  of  x  out  to  one  mile  for  the  five  stability  conditions  used  in 
Appendix  A,  the  lengths  having  been  converted  to  feet.  Appendix  B  gives  an  example 
of  the  use  of  this  graph  in  a  hypothetical  planning  problem. 
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Appendix  A 


This  Appendix  contains  graphs  of  dosage  divided  by  line  source  strength  (D/Q) 
in  minutes  per  square  foot  vs.  distance  from  the  release  line  (x)  in  feet.  For  each 
of  the  five  stability  categories  listed  in  Table  1,  graphs  of  D/Q  vs.  x  are  given  for 
assumed  release  heights  of  0,  50,  100,  150,  200,  250,  300,  400,  and  500  feet  above 
ground.  Thus  there  are  45  separate  graphs.  Each  graph  shows  4  lines  correspond¬ 
ing  to  solutions  for  transport  wind  speeds  of  5,  10,  20,  and  30  knots.  Estimates  for 
other  wind  speeds  can  be  found  by  multiplying  the  value  of  D/Q,  at  10  kt,  by  the  ratio 
10  kt/actual  wind  speed.  (For  further  details,  consult  the  test.) 
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Appendix  B 


As  an  examine  of  how  the  graphs  presented  here  might  be  used  in  planning  a 
particular  operation,  let  us  consider  the  following  problem.  Assume  that  a  grid  of 
samplers  has  been  laid  out  in,  say,  a  square  area  1  mile  on  a  side.  Also  assume 
that  a  line  source  200  feet  above  the  surface  is  to  be  emitted  in  moderately  un¬ 
stable  conditions.  It  is,  of  course,  desired  that  the  maximum  dosage  be  found  with¬ 
in  the  array  of  samplers.  The  first  problem  may  be  to  determine  how  far  from  the 
edge  of  the  grid  the  line  should  be  laid  down  so  that  the  maximum  will  be  within  the 
grid.  Assuming  the  wind  direction  is  parallel  to  one  of  the  sides  of  the  grid,  we 
would  first  check  the  appropriate  graph  in  Appendix  A.  The  graph  for  a  moderately 
unstable  condition  and  a  release  height  of  200  feet  gives  a  maximum  about  3500  feet 
from  the  release.  To  be  fairly  sure  we  have  found  the  maximum,  it  appears  that  we 
should  release  even  farther  back.  Therefore,  a  release  about  1  mile  from  the  down¬ 
wind  edge  would  be  indicated. 

The  next  problem  would  be  to  determine  the  length  of  the  release  line.  The  value 
of  x  has  already  been  set  at  1  mile.  Figure  2  indicates  that,  for  moderately  un¬ 
stable  conditions,  the  value  of  L-Lq  should  be  about  1750  feet.  Since  we  would  like 
to  ensure  that  a  significant  number  of  samplers,  at  the  farthest  distance  from  the 
source,  are  sampling  in  conditions  that  simulate  an  infinite  line  source,  a  value  of 
Lq  equal  to  1/2  mile  or  2600  feet  would  be  appropriate.  Therefore  L,  die  length  of 
the  release  line,  should  be  1750  +  2600  or  4350  feet. 

Now  a  final  check  on  the  feasibility  of  this  operation  should  be  made.  Assume 
an  aircrift  is  to  lay  out  the  particular  line  source  whose  length  should  be  4350  feet. 
Since  the  tanks  from  which  the  tracer  will  be  emitted  will  have  some  maximum 
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capacity,  say  G  grams  of  tracer,  then  the  maximum  line  source  strength,  Q',  is 
fixed  at  a  value  of  G/4350  grams/foot.  Also,  the  samplers  will  have  a  minimum 
dosage,  D',  below  which  they  will  not  give  representative  values.  Thus,  a  value 
D'/Q'  can  be  found  below  which  the  data  will  be  useless.  One  can  then  go  back  to 
the  graph  of  D/Q  vs.  x  to  see  if  values  greater  than  D'/Q'  would  be  expected  to 
cover  a  sufficient  area  of  the  grid  to  ensure  adequate  definition  of  the  distribution 
of  tracer.  If  adequate  definition  of  the  tracer  is  not  possible  then  alternatives, 
such  as  flying  lower  or  flying  more  than  once  along  the  line  while  the  samplers  are 
left  running,  must  be  considered. 
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